Phase diagram of two-color quark matter at nonzero baryon and isospin density 
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We investigate the properties of cold dense quark matter composed of two colors and two flavors 
of light quarks. In particular, we perform the first model calculation of the full phase diagram 
at nonzero baryon and isospin density, thus matching the model-independent predictions of chiral 
perturbation theory at low density to the conjectured phase structure at high density. We confirm 
the presence of the Fulde-Ferrell (FF) phase in the phase diagram and study its dependence on 
the tunable parameter in the Lagrangian that simulates the effects of the quantum axial anomaly. 
As a byproduct, we clarify the calculation of the thermodynamic potential in the presence of the 
FF pairing, which was previously based on an ad hoc subtraction of an unphysical cutoff artifact. 
Furthermore, we argue that close to the diquark (or pion) Bose-Einstein condensation transition, 
the system behaves as a dilute Bose gas so that our simple fermionic model in the mean-field 
approximation is not quantitatively adequate. We suggest that including thermal fluctuations of 
the order parameter for Bose-Einstein condensation is crucial for understanding available lattice 
data. 
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I. INTRODUCTION 

The attempts to understand the phase structure of 
quantum chromodynamics (QCD) together with the in- 
famous sign problem at nonzero baryon chemical poten- 
tial, have triggered interest in several QCD-like theories. 
These include QCD at nonzero isospin density QCD 
with adjoint quarks and two-color QCD Q. The common 
feature of these theories is that they are free of the sign 
problem. This on the one hand admits their straight- 
forward simulation using lattice Monte-Carlo techniques 
and a check of these methods against mo del- independent 
predictions of chiral perturbation theory. On the other 
hand, it allows a direct test of various model approaches 
at high temperature and/or density where chiral per- 
turbation theory is not applicable. Hence the study of 
QCD-like theories contributes to our understanding of 
the physics of strongly-coupled gauge theories at nonzero 
temperature and density. 

In this paper we investigate the thermodynamics of 
two-color QCD, which has attracted most interest among 
the QCD-like theories. We focus on the case of two light 
quark flavors with nonzero baryon and isospin chemical 
potentials. Although the underlying mechanisms of con- 
finement and chiral symmetry breaking in two-color QCD 
are believed to be the same as in the real, three-color 
QCD, the actual physical behavior is quite different, es- 
pecially at nonzero baryon density. First, two quarks can 
combine into a color singlet (leaving the gauge symme- 
try unbroken) which implies that at high density where 
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Cooper pairing near the Fermi surface occurs, the system 
behaves as a baryon superfluid rather than a color super- 
conductor. Second, even in the confined phase, the color- 
less baryons are formed by two quarks, which means that 
"nuclear" matter is actually a Bose gas and finite baryon 
density is realized through Bose-Einstein condensation 
(BEC). 

The very fact that two-color QCD is free of the sign 
problem follows from the (pseudo)reality of the gauge 
group representation in which the quarks transform. This 
has another notable consequence: the Pauli-Gursey sym- 
metry U H| connecting quarks with antiquarks. As a 
result, the spectrum of Nambu-Goldstone (NG) bosons 
of the spontaneously broken global flavor symmetry in- 
cludes pseudoscalar mesons as well as diquarks in a sin- 
gle multiplet. In fact, it is the presence of light particles 
carrying baryon number that allows to study the phase 
diagram of two-color QCD at nonzero temperature and 
baryon density using the techniques of effective field the- 

ory nag. 

While chiral perturbation theory provides a model- 
independent parameterization of all observables in terms 
of a few low-energy couplings, its validity is limited to 
the region of low density. Several approaches have there- 
fore been used to extend its results. These include 
the linear sigma model [7], random matrix theory 
or high-density effective theory The Nambu-Jona- 
Lasinio (NJL) model was first applied to two-color QCD 
in Ref. [1(| for one quark flavor and in Ref. 11 1 for two 
flavors. In [l2[ the crossover from BEC of molecular di- 
quarks to the Bardeen-Cooper-Schrieffer (BCS) pairing 
of quarks was studied. Nonzero isospin chemical poten- 
tial imposes stress on this pairing and when sufficiently 
large, induces a first-order transition to the normal phase. 
Around this transition, a narrow window with inhomo- 
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geneous pairing may exist, as conjectured in Q and ex- 
plicitly verified in jl3j . Despite the number of papers on 
the subject pursuing different analytic approaches, the 
main body of work has been done on the lattice, using 
both the strong-coupling expansion [1J] and Monte-Carlo 
simulations [TBI. fl6l] . 

The plan of the paper is as follows. In Sec. [TT] we in- 
troduce the NJL model for two-color QCD and demon- 
strate how it can be derived in a way which makes the 
extended symmetry of two-color QCD manifest. Section 
IIIII is entirely devoted to the mean-field approximation 
within this model. The reason is that inclusion of the 
Fulde-Ferrell (FF) pairing brings nontrivial regulariza- 
tion issues, in particular in connection with the widely 
used sharp three-dimensional cutoff. Section ITVl presents 
the numerical results for the phase diagram at nonzero 
baryon and isospin density. We include chiral, diquark 
and pion condensates with the possibility of spatial (FF) 
modulation of the latter two. Using a generalized NJL 
interaction, we also allow for isospin breaking by dif- 
ferent chiral condensates of the u and d quarks. We 
confirm the presence of the FF phase at high chemical 
potential, in agreement with We show that the 

region in the phase diagram where FF pairing occurs 
depends sensitively on the parameter in the Lagrangian 
that measures axial symmetry breaking by the quantum 
anomaly. In Sec. [V] we investigate the collective modes in 
the NJL model. We calculate their propagators and de- 
termine their masses and dissociation temperatures. In 
Sec. lVII we discuss in some depth two-color quark matter 
at nonzero baryon and zero isospin density, which corre- 
sponds to the conditions being simulated on the lattice. 
Recent progress in the lattice computations led to first 
attempts to probe BEC of diquarks and the region of 
moderate baryon density One of the motivations 

for the present work was an attempt to understand the 
available lattice data on the thermodynamics of cold and 
dense two-color quark matter. Even though we have 
not been able to demonstrate quantitative agreement be- 
tween the lattice simulation and model calculations, we 
outline a possible direction along which it can presum- 
ably be reached. Finally, in Sec. IVIIl we summarize and 
conclude. 



II. NJL LAGRANGIANS FOR TWO-COLOR 
QUARK MATTER 

In any model attempting to mimic some features of full 
QCD, it is mandatory that its Lagrangian has the same 
symmetries as the underlying theory. In 11| the NJL in- 
teraction for two-color QCD was derived by a Fierz trans- 
formation of the four-quark interaction inspired by one- 
gluon exchange, in which the flavor symmetry is manifest. 
Here we introduce a formalism which makes it possible to 
construct NJL interactions with desired symmetry prop- 
erties directly. Our discussion will parallel that of Kogut 
et al. Yet we believe it is worthwhile to show some 



details, since instead of working with the Weyl spinors for 
the quark fields as in Q, we keep the Dirac spinors and 
introduce the Nambu doublet notation, as is common in 
the literature on the NJL model. 

Our starting point is the Lagrangian for the quark sec- 
tor of two-color QCD, 



2cQCD 



(1) 



where D^i/j = (<9 M — iga a A^)^p is the color-SU(2) co- 
variant derivative, mo is the current mass common to 
all quark flavors, and the flavor indices are suppressed. 
Next we replace the usual Dirac spinor, consisting of 
the left- and right-handed components, with the purely 
left-handed Nambu spinor and its (conveniently defined) 
charge conjugate, 



* = 



-,-rT 



where ip c — C^p' is usual Dirac charge conjugation. The 
two-color QCD Lagrangian ([T]) can then be rewritten in 
the block-matrix form 



2cQCD 
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* + H.c. 



(2) 



All entries of the 2x2 Nambu mass matrix are themselves 
matrices in flavor space. Without going into details which 
are covered in literature [2j, we note that thanks to the 
(pseudo)reality of the gauge-SU(2) generators, J?2cQCD 
with Nf quark flavors has an apparent SU(27V/) sym- 
metry in the chiral limit. The baryon number, being 
generated by 



R 1 { 1 

B - 2 l0 -1 



is incorporated as a part of this simple group, while a 
mere change of the phase of VE' corresponds to axial U(1)a 
transformations. This symmetry is known to be broken 
at the quantum level by the axial anomaly. 

As is clear from Eq. (JSJ), the mass operator in the 
Nambu notation is proportional to 







The standard chiral condensate which represents the 
ground state in the vacuum has the same flavor structure 
as the mass term. Therefore, the multiplet of mesons, one 
of which is to condense in the ground state and the oth- 
ers to become NG modes, transform as a rank-2 antisym- 
metric tensor of SU(27V/). For Nf = 2 we make use of 
the local group isomorphism SU(4) ~ SO (6) to conclude 
that with a suitable basis of 4 x 4 antisymmetric matrices, 
£ = {Sj}| =1 , the operator ^E ,C E^E' transforms as a com- 
plex vector of SO (6). The basis matrices are unitary and 
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satisfy the orthogonality condition X^E 
They may be chosen to carry the quantum numbers of 
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the individual meson modes: three pions, sigma, diquark 
and the antidiquark. The explicit form of the matrices 
with these properties was given in fTrl j . 

It is now obvious how to construct a NJL interaction 
Lagrangian with SO(6) symmetry. As long as we only 
consider the Lorentz-scalar bilinear operator \E ,C S5', thus 
selecting bosonic modes with specific quantum numbers 
[l8j . there are only two ways to make an SO(6) invariant 
out of the vector *f rC £\E r : taking a scalar product with 
itself or its complex conjugate, 



'U(4) 
-^SU(4) 



G 



(¥ c l!tf) S 



H.c. 



(3) 



The subscripts emphasize the important fact that while 
Jz?u(4) i s invariant under the U(1)a group of phase trans- 
formations of 5', Jz?su(4) obviously breaks it. Inserting 
explicit expressions for the matrices, these two inter- 
actions can be readily rewritten in terms of the original 
Dirac spinor ip. It is convenient to work with the follow- 
ing combinations, 



— 2 [-^U(4) + -2?SU(4)] 

= G (^l/>) 2 + (#y 5 T^) 2 + \Wl5(T2T 2 1p\' 
-^2 = 2 [-^U(4) — -^SU(4)] 



where f are Pauli matrices in isospin space. The quantum 
numbers of the various channels are now apparent: 
contains the sigma, pions and the scalar (anti)diquark, 
while Jz?2 contains their parity partners. 

Note that ifi is essentially the minimal two-flavor NJL 
Lagrangian [l9j , augmented by a diquark term. The only 
difference with Lagrangians used in model investigations 
of dense (three-color) quark matter is that here the SU(4) 
symmetry fixes the coupling in the diquark channel to be 
the same as for the meson terms. This is the Lagrangian 
used in [lll[l~2j. However, since we have two independent 
interaction terms in Eq. ([3]) , the most general interaction 
one can construct is a linear combination 



-Spirit — -2*1 + C-^2j 



(4) 



and was written down for the first time in [20( . At £ = 1 
the Lagrangian reduces to ifu(4)i preserving the axial 
U(1)a, and is a direct descendant of the current-current 
interaction used in 11] . The case of ( — —1, that is, 
Jzfsu(4), corresponds to the two-flavor instanton-induced 
interaction, breaking U(1)a maximally |4j. In the present 
paper £ is treated as a free parameter, the only restriction 
being £ < 1. This follows from the fact that in the vac- 
uum, the scalar chiral condensate must be preferred over 
the pseudoscalar one, as guaranteed by the Vafa-Witten 
theorem [2l[. We note in passing that in the three- 
color NJL model, one usually introduces the parameter 



a = (l — C)/2 which measures explicit U(1)a breaking by 
the axial anomaly so that the interaction Lagrangian can 
also be expressed as Jzf; n t = (1 — a)Jz?u(4) + a -^su(4)- The 
typical values of a lie somewhere in the range 0.1 — 0.2 

The full NJL Lagrangian now reads 

if = 1p[i0 - TO + \lf){p.B + T 3 Hi)]lp + JS?int) 

where ps and /i/ are baryon number and isospin chemi- 
cal potentials, both being assumed positive without lack 
of generality. In 20] this Lagrangian with the generalized 
interaction (jj) was employed to study the masses of the 
scalar mesons and hence the axial symmetry restoration 
in the meson spectrum. Here it will be crucial at high 
baryon and isospin chemical potentials, for it allows to in- 
troduce different it-quark and d-quark chiral condensates 
within the mean-field approximation. 



III. MEAN-FIELD APPROXIMATION 

As mentioned above, the scalar chiral and diquark con- 
densates are preferred to the pseudoscalar ones. Since 
we will not be interested in the corresponding non-NG 
modes (eta meson and the pseudoscalar diquark) either, 
we will completely neglect the interaction terms in the 
t/jij^ip and ?/; c <7272'0 channels. The remaining four chan- 
nels can be rewritten using the Hubbard-Stratonovich 
transformation. The Lagrangian can then be written as 



if = tp[i0-mo + ^ o (fi B +T3^i)-(T-i^ 5 T-TT-Cf-p] ip 

(5) 



+ 5(AVtswh^Hx.) - ^(<j 2 +7T 2 + (p 2 + |A| 2 ), 



where a, if, p, A are the collective fields with the quan- 
tum numbers of a scalar isoscalar, pseudoscalar isovector, 
scalar isovector, and scalar diquark, respectively. 

In the mean-field approximation the collective fields 
are replaced with their vacuum expectation values. The 
presence of both baryon number and isospin chemical po- 
tentials induces a mismatch between the Fermi levels of 
quarks of different flavors. When the mismatch is too 
large, it breaks the BCS-like pairing and the system goes 
over to the normal phase via a first-order phase transi- 
tion. Around this transition, a narrow window may ex- 
ist where pairing with nonzero total momentum occurs. 
Such kind of pairing, well known in condensed matter 
physics, was first considered in the context of dense quark 
matter in [23| . see also [24j for a comprehensive review. 

In literature on color superconductivity of dense 
(three-color) quark matter (see [25[ for recent reviews), 
there are in fact several other candidates for the state 
of dense matter in presence of a Fermi surface mis- 
match that compete with the inhomogeneous (FF) pair- 
ing. These includephases with a gluon condensate [26j j . 
meson condensate [27| , or a meson supercurrent [28j | . Let 
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us explain briefly why we do not consider these possibil- 
ities in the present context. First, in a two- flavor (three- 
color) color superconductor, the single plane wave FF 
state is equivalent to the BCS state with a uniform gluon 
condensate; they are connected by a gauge transforma- 
tion (26J. However, this is not the case in two-color QCD 
where the diquark pairs are color singlets. For the same 
reason, the gluon condensate will be disfavored since it 
induces nonzero color charge (which in three-color QCD 
compensates for the charge of the Cooper pairs) . 

The possibility of a meson condensate strongly de- 
pends on the meson spectrum, and hence on the pattern 
of global symmetry breaking. In three-color QCD me- 
son condensation is likely to occur in three-flavor matter 
where there is an octet of light pseudoscalar mesons [13] ■ 
On the other hand, in two-color QCD with two flavors, 
considered in this paper, the only mesons present in the 
low-energy spectrum of the diquark condensation phase 
are the pions. It can be shown that at pi = 0, the pion 
mass is equal to ps [HI- This eventually does lead to 
pion condensation, but only when pi > ps- There is 
no meson to condense inside the diquark condensation 
phase, thus leaving the FF pairing the only candidate for 
resolving the Fermi surface mismatch. 

In the present paper, we will therefore take into ac- 
count the possibility of inhomogeneous diquark and pion 
condensates. To that end, note that the Lagrangian of 
two-color QCD ([TJ is invariant (among others) under 
the discrete symmetry of charge conjugation of <i-quarks. 
The reason why we can charge conjugate one component 
of the quark flavor doublet and still keep the Lagrangian 
gauge invariant is once again the (pseudo)reality of the 
doublet representation of the gauge-SU(2). This discrete 
symmetry remains intact even in presence of the baryon 
number and isospin chemical potentials, provided they 
are interchanged simultaneously with the charge conju- 
gation. As a consequence, the phase diagram of two-color 
QCD in the (pb,Pi) plane is symmetric with respect 
to interchanging of the axes accompanied by the simul- 
taneous interchanging of the diquark and pion conden- 
sates Moreover, the diquark and pion condensates 
never occur simultaneously, the former being favored for 
fJ-B > fJ>i, while the latter for ps < pi- This is rather 
obvious at high chemical potentials where quarks form a 
Fermi sea: at ps > fJ-i a Fermi sea of <i-quarks exists, 
thus making formation of ud pairs, and hence diquark 
condensation, possible. Likewise, at pb < Hi, (i-quarks 
form a Fermi sea, leading to ud Cooper pairs and pion 
condensation. At low chemical potentials, the mutual 
exclusion of the pion and diquark condensates may be 
verified using chiral perturbation theory. 

With the above argument in mind, we will set 7? = 
and consider the region < pi < ps- The full phase 
diagram is easily constructed by reflection with respect 
to the line ps = fJ-i, accompanied by the replacement of 
the diquark condensate with the pion one. For the coor- 
dinate dependence of the diquark condensate we will as- 
sume the simplest, single plane wave form, Ae 2tq x with 



real A. This is the FF pairing [30j]. As a side remark 
let us add that the FF pairing is not really the energet- 
ically most favored state. The energy may be further 
decreased by forming a more complicated spatial struc- 
ture such as the truly inhomogeneous standing wave-like 
Larkin-Ovchinnikov state 3l| , or even various crystalline 
structures [32]. When the periodic structure of the or- 
der parameter is not limited to a linear combination of a 
finite number of plane waves, the first order transition be- 
tween the BCS pairing and the FF phase is weakened and 
eventually becomes second order, enlarging the region in 
the phase diagram occupied by the inhomogeneous phase 
[33( | . FF pairing in the context of pion condensation was 
investigated before in 34 1. 

All in all, we take into account the following mean 
fields: a, p (standing for the third component of p), and 
/S.e 2lq x . Note that p gives rise to splitting of the con- 
stituent quark masses, being M u ^ = ttiq + er ± (^p. The 
coordinate dependence of the diquark condensate may 
be removed by a unitary transformation at the price of 
modifying the kinetic term of quarks. After integrating 
out the quarks and performing the Matsubara summa- 
tion, the thermodynamic potential density of the model 
becomes 

fi(cr, p, A, o) 1,2 j. 1 . o, 

V -4G^ +C f +A) (6) 

^Trlog(l + e-^). 



(2tt) 3 



Here J#k is the 16 x 16 matrix Hamiltonian whose block 
elements in the Nambu space (i/} r ,(T2tpq), read 



P4)n =a-(k + q)+ 7o (M + ( P t 3 ) 

(j#k)l2 = - 7075A, 
(^fc)21 = +7075A, 

{M)22 =a-(k-q) + 7o(M - Cpr 3 ) 



:(Pb +r 3 p,i), 



where we introduced the notation M — mo + a for the 
average constituent quark mass, and a — 707 are the 
usual Dirac matrices. Note that the zero-point energy 
contribution is missing in Eq. (j6|) since TrJfk = 0. The 
eigenvalues of the Hamiltonian, representing the energies 
of fermionic quasiparticles, cannot be found analytically. 
When calculating the thermodynamic potential, one has 
to resort to numerical diagonalization of the Hamilto- 
nian. In [l3| this technical difficulty was avoided by using 
an approximate analytic expression for the quasiparticle 
dispersion relations. In the present paper the thermo- 
dynamic potential is evaluated without introducing any 
further approximation. 

The naive expression for the thermodynamic potential 
([6]) is actually ill-defined because the momentum integral 
is badly divergent. The integral will be regulated with 
a three-momentum cutoff, A, as usual. However, with 
nonzero FF momentum q one must be very careful at 
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what stage of the calculation the cutoff is applied. Try- 
ing to impose it naively to Eq. ^ leads to unphysical 
artifacts including a q-dependent quadratic divergence. 
This was previously removed by an ad hoc subtraction 
[IH, [35j , making sure that the thermodynamic potential 
is independent of q once A = 0, which follows from the 
very definition of the FF order parameter. We would 
like to argue here that such a subtraction follows natu- 
rally when the thermodynamic potential is first properly 
regularized, and only after then the unitary transforma- 
tion leading to Eq. © is performed. Since the argu- 
ment is rather general and applies in particular also to 
the analogous nonrelativistic problem, we present it sep- 
arately in the Appendix. The final expression for the 
subtracted thermodynamic potential, using the function 
fl((X, p, A, q) defined by Eq. ([5]), is 

O sub (o-, p, A, q) =ti(a, p, A, q) 

The addition and subtraction operations are understood 
to be performed inside the integration. This prescription, 
differing slightly from that used in [l3[ , also ensures that 
the thermodynamic potential is independent of q once 
A = 0. Moreover, the thermodynamic potential Sl su b 
has a well defined minimum which determines the equi- 
librium state. 



IV. PHASE DIAGRAM AT NONZERO BARYON 
AND ISOSPIN DENSITY 

The phase diagram of two-color QCD with two quark 
flavors in the (un, pi) plane and at zero temperature was 
first drawn in [5J . At low chemical potential its structure 
can be determined in a model-independent manner using 
the chiral perturbation theory. One thus finds that when 
both chemical potentials are smaller than m^, the vac- 
uum mass of the pion-diquark multiplet of pseudo-NG 
bosons of the spontaneously broken SU(4) flavor symme- 
try, the system is in the vacuum phase. Only the chiral 
condensate is nonzero, equal to its vacuum value inde- 
pendently of the chemical potential. As soon as one of 
the chemical potentials, say pb (remember the reflec- 
tion symmetry of the phase diagram), exceeds m T , the 
mode carrying the associated charge — in this case the 
diquark — condenses via a second-order phase transition. 
Since this picture relies just on the symmetry and the 
pattern of its spontaneous breaking, it is not surprising 
that results obtained within chiral perturbation theory, 
lattice simulations [l5[, and model calculations [TTJ] are 
in good quantitative agreement. At moderate chemical 
potentials one therefore finds BEC phases with a diquark 
(for pb > pi) or pion condensate (for ps < pi), sepa- 
rated by a first-order phase transition line at pb = pi- 

On the other hand, at high chemical potential one ex- 
pects the physics of the system to be dominated by the 
Fermi sea of quarks. Spontaneous symmetry breaking 



is realized by a formation of Cooper pairs via the stan- 
dard BCS scenario. Since the Cooper pairs carry the 
same quantum numbers as the condensing modes in the 
low-density BEC phase, there is no phase transition be- 
tween the two regimes despite the fact that they feature 
qualitatively different physical behavior. One speaks of 
the BCS-BEC crossover |3q[. I n the context of two-color 
QCD, this was studied in [12j; we will comment on it 
later in Sections IVl and IVT1 

The qualitative shape of the phase diagram at high 
chemical potential was conjectured in [5j based on the 
experience with real, three-color QCD at weak coupling. 
The essence of the argument is as follows. Close to the 
Pb and pi axes, one will find diquark and pion conden- 
sates as predicted by BCS theory. However, since both of 
them involve pairing of quarks of different flavors, they 
will become less and less favored towards the "diagonal" 
of the phase diagram; pi generates a mismatch between 
the Fermi levels of u and d quarks in the diquark conden- 
sation phase, while ps plays the same role in the pion 
condensation phase. When the mismatch becomes large 
enough, the BCS pairing is no longer energetically favor- 
able and it breaks in a first-order phase transition. In 
a narrow window around this transition, a phase where 
Cooper pairs carry nonzero total momentum may exist. 
Even after the cross-flavor pairing is made impossible, 
quarks of the same flavor can still pair, since there is 
no Fermi level mismatch by definition. Such pairs are re- 
quired to carry nonzero spin or orbital momentum by the 
Pauli principle 10]. In a narrow strip along the diagonal 
of the phase diagram, the chemical potential of d-quarks 
is almost zero, and they should therefore be confined and 
exhibit the one-flavor chiral condensate. Thus, only the 
it-quarks will pair in this region. 

We employed the NJL model to verify the conjecture 
made in [5j . In order to obtain concrete numbers, we used 
the two-color NJL parameter set, established in [20] : G = 
7.23 GeV~ 2 , m = 5.4 MeV, and A = 657 MeV. This 
corresponds to the physical pion mass = 140 MeV, 
pion decay constant /„■ = 75.4 MeV, and (one-flavor) chi- 
ral condensate (^> u ip u ) = (-218 MeV) 3 . The latter two 
quantities were obtained from usual three-color values by 
rescaling with respect to the number of colors. The pa- 
rameter £ is treated as free in this paper. 

While we included the FF-type diquark and pion con- 
densates as well as independent chiral condensates for u 
and d quarks, we did not, for the sake of simplicity, take 
into account spin-one pairing. This would most likely 
further split the region in the upper right corner of the 
phase diagram into domains with or without spin-one 
pairing of individual flavors. On the other hand, it is not 
expected to alter the structure that we find. The rea- 
son is that spin-one pairing is usually characterized by a 
small pairing gap as compared to spin-zero pairing, and 
the pairing moreover does not occur isotropically on the 
whole Fermi surface. Consequently, the energy gain from 
spin-one pairing is minute so that it cannot compete with 
phases with other types of pairing. 
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FIG. 1. Phase diagram of two-color QCD with two quark 
flavors in the (pb,Pi) plane at zero temperature. Various 
phases are labeled by the corresponding nonzero condensates. 
The p condensate in the upper right corner only appears at 
( / 0. Thick and thin lines denote first- and second-order 
transitions, respectively. The black solid lines were obtained 
with f = (a = 0.5), the red dash-dotted lines with £ = 0.8 
(a — O.f), and the blue dashed lines with £ = 1 (a = 0). The 
position of the transitions at high chemical potentials only 
changes appreciably for £ close to one. The £ = 0.6 lines 
would almost coincide with the £ = ones shown. 



The phase diagram at zero temperature obtained with 
the above set of parameters and condensates taken into 
account is shown in Fig.[TJ One can see that at moderate 
chemical potentials, there is no FF phase and the system 
goes directly to the phase with just chiral condensate(s). 
We avoid calling this a normal phase, because the size of 
the chiral condensate(s) depends sensitively on the pa- 
rameter £. For £ = 0, p = and the common chiral 
condensate of the u and d quarks is indeed very small. 
As ( increases, the flavor chiral condensates (J u ,d — &±Cp 
split, or in other words, p becomes nonzero. Naturally, 
a u remains very small since the u-quark chemical poten- 
tial is very large, while ad grows towards the diagonal 
of the phase diagram. In fact, in the U(4) symmetric 
limit (£ = 1) the mean- field gap equations for a u and ad 
in absence of diquark and pion condensates completely 
decouple. The d-quark condensate a^ is then at the di- 
agonal of the phase diagram (and around) equal to its 
vacuum value. As soon as the U(1)a symmetry is ex- 
plicitly broken, the flavors become entangled and the d- 
quark condensate drops fast with increasing a. However, 
despite the strong dependence of the flavor chiral conden- 
sates on C, the position of the phase transition does not 
change much. Some shift is actually only visible as C, ap- 
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FIG. 2. The high chemical potential corner of the phase di- 
agram. The notation for the lines is the same as in Fig. \T\ 
The region where the momentum q is nonzero is indicated as 
"FF" . The additional thick blue dashed line is associated with 
the first-order chiral restoration transition at £ = 1, explained 
in the text. Note that the chemical potential used to label the 
horizontal axis is the baryon number one so that the corre- 
sponding quark chemical potential in the range displayed is 
still below the cutoff A. Hence, cutoff artifacts, though non- 
negligible, will not alter the qualitative structure of the phase 
diagram. 



proachcs one. Note that for the sake of clarity, we do not 
indicate the BCS-BEC crossover in Fig. [TJ This appears 
at pb,i ~ 229 MeV; its position is to a very good pre- 
cision determined by the expression (2Moto^) 1 / 3 , where 
Mo is the constituent quark mass in the vacuum (T^ |. 

The domain of high chemical potentials is shown in 
detail in Fig. [2j The FF phase appears only at rather 
high chemical potentials. This is in agreement with the 
observation that at in strongly coupled systems the re- 
gion with the FF phase is diminished [H, H3| • In quark 
matter one cannot tune the coupling arbitrarily. Never- 
theless, the weakly coupled BCS pairing is achieved ef- 
fectively at high chemical potential since the ratio A / ps 
drops with increasing chemical potential. Even though 
we used a different parameter set than in [13] and a more 
refined method to evaluate the thermodynamic potential 
at nonzero FF momentum q, our results are in reasonable 
quantitative agreement. 

As £ increases, the FF phase is obviously expelled from 
the phase diagram. This is a direct consequence of the 
fact that the energy gain from FF pairing is rather small 
due to both the small size of the pairing gap and the 
fact that pairing occurs only on a small part of the Fermi 
surface. It therefore becomes energetically more favor- 
able to form a state with a large splitting of constituent 
u and d quark masses, which in turn inhibits the cross- 
flavor pairing. On the other hand, one should note that 
in the FF phase, p does take a nonzero, albeit small 
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FIG. 3. Meson masses as a function of temperature at zero 
chemical potential. Solid line: common mass of pions and 
the (anti)diquark. Dash-dotted line: sigma mass. The thin 
dashed line indicates the position of the two-particle contin- 
uum, that is, twice the constituent quark mass. 



value. The reason is that the difference of the constituent 
quark masses partially compensates for the chemical po- 
tential mismatch, thus reducing the mismatch of Fermi 
momenta. 

Fig. [2] also features an additional first-order transition 
line which occurs only for £ = 1 and is defined by the 
condition \Ib — = const. This can be explained using 
the observation made above that at ( — 1 the gap equa- 
tions for a u and ad completely decouple. The two flavors 
then live their separate lives, and instead of (J,b,i it would 
be more convenient to use the basis fi u ,d — \{l^B ± /ij). 
The first-order line in question corresponds to the critical 
value of fid at which a strong first-order chiral symme- 
try breaking/restoration phase transition occurs. The 
critical value, fid ~ 330 MeV, extracted from the phase 
diagram in Fig. [2] , is in excellent agreement with an el- 
ementary calculation within a one-flavor version of our 
model (O, keeping the chiral condensate ad as the sole 
order parameter. 

Once £ drops below one, that is, axial U(1)a symmetry 



is explicitly broken by nonzero a, the two quark flavors 
become entangled. This smears the phase transition in 
the d-quark sector, making it a smooth crossover. We 
have checked numerically that there is no sharp phase 
transition already for very small a, and for the case 
a = 0.1 displayed in Fig. [5] by red dash-dotted lines, 
the crossover is already rather gradual. 



V. COLLECTIVE MODES 

With the semi-bosonized Lagrangian ([5]), the propaga- 
tors of the bosonic collective modes are straightforwardly 
derived by integrating out the quarks and taking the sec- 
ond variational derivative with respect to the collective 
fields Xi- The general expression for the inverse propa- 
gator at imaginary (Matsubara) frequency i£In is [2C| 



x Tr 



d 3 k 

(27)3 



dxi i(w„ + Oat) - d\j iu n - Jl? k „ 



This formula is valid for the a, 7? , A modes we are inter- 
ested in. For modes from the other chiral multiplet, a 
factor C must be added to the constant part of the prop- 
agator. In general, when some of the condensates break 
the symmetry spontaneously, different modes may mix. 
We will only give explicit expressions for the propagators 
in the case that just the isospin-singlet chiral conden- 
sate (a) is nonzero. This covers the connected domain in 
the phase diagram containing the origin. Since the BEC 
transition to either diquark or pion condensed phase is of 
second order, the knowledge of these propagators allows 
one to locate these transitions by means of the Thouless 
criterion [38| . For the sake of simplicity of the result- 
ing formulas, we will set the external three-momentum 
to zero. (A more complete analysis of collective modes 
may be found in [20| ; the evolution of the collective mode 
spectrum as a function of baryon chemical potential has 
also been studied on the lattice |39(.) With these sim- 
plifications in mind, the inverse propagators of the six 
modes of interest to us read, after analytic continuation 
to general complex frequency u>, 
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d 3 k i-mp-Mv) 

(2ir) 3 ui + fij + 2eek 
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(8) 



We use the usual notation for the Lorentz-invariant dis- 
persion relation, = \/k 2 + M 2 , and the Fermi-Dirac 
distribution function, f(x) = \ j(eP x + 1). Furthermore, 
= e fc + f (a*-B + A*/) denotes the in-medium energy of 
the u-(anti)quark, and Cfe = e fe + §(MB — Ml) that of the 
d-(anti)quark. The propagators of tt~ and of the antidi- 
quark are obtained from those of n + and the diquark by 
changing the sign of [ii and fin, respectively. 

In Fig.[3]we plot the meson masses as a function of tem- 
perature at zero chemical potential. Note that thanks to 
the unbroken SO(5) symmetry, the pions and diquarks 
are exactly degenerate even away from the chiral limit. 
The masses are determined from the zeroes of the real 
part of the inverse propagators ([8]) . As long as the mode 
is stable, this definition coincides with the physical, pole 
mass. Once the two-particle threshold is crossed, the defi- 
nition of the mass becomes ambiguous and it is more con- 
venient to discuss the full spectral density, p(uj). This is 
obtained from the imaginary part of the retarded Green's 
function, and is related to the inverse propagators © by 

_ ImD' 1 ^ + ie) 

P(UJ ' ~ ~ [RcD-^u + ie)} 2 + [ImD^iu + ie)} 2 ' 

The spectral densities as a function of fre- 
quency for five different values of temperature 
(140, 180, 220, 260, 300 MeV) are plotted in Fig. H One 
can see clearly that all modes constitute well-defined 
quasiparticles below the chiral restoration crossover; the 
pions and diquarks are bound and therefore stable, while 
the sigma mass is just above the threshold and thus 
gives rise to a very sharp peak in the spectrum. Once 
the temperature increases, the constituent quark mass 
rapidly drops, the level of the two-particle continuum 
consequently decreases, and all mesons become broad 
and eventually more or less washed out of the spectrum. 
For temperatures above about 300 MeV the notion of 
mass becomes meaningless. 

The range of stability of the modes is more precisely 
determined by the so-called dissociation temperature, de- 
fined as the temperature at which the quasiparticle pole 
in the propagator (|SJ) enters the two-particle continuum 
[40j . (Note that a different definition of the dissociation 
temperature was used in [Hj].) In Fig. [5] the phase di- 
agram at zero isospin chemical potential is drawn with 
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FIG. 4. Spectral densities of the pions and diquarks (up- 
per panel) and sigma (lower panel) as a function of fre- 
quency at zero chemical potential. The various lines corre- 
spond to temperature of: 140 MeV (black, solid), 180 MeV 
(red, solid), 220 MeV (green, dashed), 260 MeV (blue, dash- 
dotted), 300 MeV (magenta, dotted). 



the dissociation temperatures of the pions, the diquark 
and antidiquark indicated. We can see that our definition 
of the dissociation temperature nicely coincides with the 
position of the BCS-BEC crossover, to be discussed in 
more detail in the following section. It is suitable to re- 
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FIG. 5. Phase diagram of two-color QCD at zero isospin 
chemical potential. The thick black solid line denotes the 
second-order diquark condensation transition. The thin black 
dashed line indicates the position of the BCS-BEC crossover, 
denned by the condition M = /is/2. The dissociation tem- 
peratures of the collective modes are represented by the: red 
dashed line for the diquark, green dash-dotted line for the 
antidiquark, blue dotted line for the pions. The region where 
the physics is dominated by the chiral condensate is labeled 
as X SB. 



mark here that when the effects of confinement are taken 
into account, for example by coupling the NJL model to 
the Polyakov loop, the transition temperature becomes 
somewhat higher than in Fig. [5] [20J. However, the qual- 
itative structure of the phase diagram is not changed. 
Also, it should be stressed that adding the Polyakov loop 
would only modify the results of this section. It would 
have no effect at all on the phase diagram in the (fiB, fJ-i) 
plane at zero temperature, discussed in Sec. IIV1 and it 
plays an absolutely negligible role in the thermodynamics 
of diquark BEC, to be investigated in Sec. I VII [16]. 
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[MeV] 

FIG. 6. Pressure and baryon density (upper panel) and energy 
density (lower panel), normalized to the Stefan-Boltzmann 
values, at zero isospin chemical potential as a function of 
baryon number chemical potential. The thick lines correspond 
to T = MeV with the following notation: the red solid line 
for pressure, the green dash-dotted line for baryon density, 
and the blue solid line for energy density. The thin lines with 
the same line notation show the results for T = 50 MeV. The 
thin black lines, almost coinciding with the zero-temperature 
NJL values, are the results of chiral perturbation theory at 
leading order. 



VI. THERMODYNAMICS OF DIQUARK 
BOSE— EINSTEIN CONDENSATION 

The main virtue of two-color QCD, providing a strong 
motivation for its study, is the possibility of lattice Monte 
Carlo simulations. It is therefore of utmost importance 
to try to match model calculations to available lattice 
data. To this end, one should note that in order to have 
a positive determinant of the Dirac operator that in turn 
makes the importance sampling possible, one needs an 
even number of quark flavors with the same chemical 
potentials. The phase diagram in the (/is,/i/) plane in- 
vestigated in Sec. IIVI is therefore not accessible to cur- 
rent lattice techniques. The lattice literature on two- 
color QCD at nonzero density invariably addresses the 
case of zero isospin chemical potential, where the indi- 



vidual quark flavors indeed share the common chemical 
potential, /is . The latest simulations [l6| study the ther- 
modynamic properties of two-color QCD with two quark 
flavors at low temperature as a function of the baryon 
chemical potential. (They also calculate the gluon prop- 
agator which, however, cannot be determined in model 
approaches based on the flavor symmetry.) 

Within our NJL model, we calculate the pressure 
from the thermodynamic potential ([7]), making an ap- 
propriate constant subtraction to ensure that the vac- 
uum pressure is zero. The baryon number density and 
energy density are in turn obtained from the pressure 
using the thermodynamic identities, tib — dp/dfiB and 
e = —p + [Ib^b + Ts, where the entropy density reads 
s = dp/dT. Since at high chemical potentials the system 
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is expected to behave as a weakly interacting Fermi gas, 
it is convenient to normalize these thermodynamic quan- 
tities by the Stefan-Boltzmann (SB) values, correspond- 
ing to a free massless Fermi gas. Taking into account 
that two-color QCD with two quark flavors includes four 
different species of Dirac fermions, these read 
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and £sb = 3psB- The numerical results for two differ- 
ent temperatures are shown by the colored lines in Fig. [6] 
TlJ- Remarkably, even for the highest chemical poten- 
tials displayed, the pressure and density are still rather 
far from the SB limit. In the mean- field approximation, 
this is obviously due to the pairing gap. 

Since near the diquark BEC phase transition, the the- 
ory behaves like a dilute Bose gas rather than a system 
of noninteracting quasiparticles, described by the NJL 
model in the mean-field approximation, we calculate for 
comparison the same quantities in a model with elemen- 
tary scalar degrees of freedom. For the sake of simplicity, 
we use the 0(6) linear sigma model at tree level (corre- 
sponding to zero temperature). This is defined by the 
Euclidean Lagrangian 
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Ha, 



(9) 



Note that the results of chiral perturbation theory Q 
are recovered by taking the limit m a oo. Also, the ex- 
pressions for the chiral condensate, a = f^m 2 /^, and 
the baryon density, ns = /is A 2 , are actually identical 
to those in the chiral perturbation theory. (The baryon 
density, nevertheless, takes on a different value due to 
the modified value of the diquark condensate.) It is in- 
teresting to observe that as the chemical potential ap- 
proaches from the right, the predictions of the linear 
sigma model and chiral perturbation theory are related 
by A 1(Tm /A xPT = m a /y/m 2 — m 2 . This means that as 
soon as the diquark condensate is nonzero, a finite sigma 
mass drives the results off the predictions of chiral per- 
turbation theory even in the limit of zero density. Given 
the above solution of the gap equations, it is straightfor- 
ward to calculate the thermodynamic quantities. Using 
again the normalization with respect to the SB values, 
one finds 



= 2Att 2 ( ^ 

PSB Kniyr 



1 



X. 



1 



1 



"lam _ 12?H / J^A f _ J_ 2 ^ 

x 2 \ mjr I \ x 4 t 



nsB 

^SB 



X 2 T 2 — 1 
2 -l 



(10) 



1 'i — 

X 2 



3x 2 



1 



<f> = (a, tt, Re A, Im A) r , where 7? is the pion triplet and 
A the complex diquark field. The parameter H breaks 
the 0(6) symmetry explicitly down to 0(5). The baryon 
chemical potential /is is incorporated by the replacement 
of the ordinary derivatives with covariant ones, acting 
nontrivially only on the diquark, DqA = (8q — /is) A, 
andA)A+ = (<9o + Ais)At. 

The tree-level analysis of this model is straightforward, 
so we proceed directly to the results. The negative mass 
squared triggers condensation of the field, which, in the 
vacuum, is aligned with the explicit symmetry break- 
ing term, er = f w = (M 2 + m 2 )/X, where the pion 
mass is determined implicitly by H — /^to 2 . As soon 
as the baryon chemical potential exceeds the pion mass, 
diquarks condense. The chiral and diquark condensates 
are then in general parameterized by a — v cos 9 and 
A = i>sin#. For the magnitude v and the angle 9 one 
finds the expressions 
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/4 
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The unphysical Lagrangian parameter M 2 may be elim- 
inated in favor of the mass of the "radial" excitation, 
conventionally identified with the sigma particle, m 2 = 
2M 2 + 3m 2 . All quantities of interest can then be ex- 
pressed in terms of observables. For instance, the diquark 
condensate is given by 



where we used the shorthand notation x = fMs/Tn^ and 
r = Too- /m-n- Similar expressions were obtained for pion 
condensation in [421 ] . It is amusing, though perhaps not 
very physical, that all the expressions (jlOp have a nonzero 
limit as x — > oo. This takes even the same value for all 
three observables, 247t 2 / 2 /(to 2 — m 2 ). 

Taking the values of m n and f v over from the NJL 
model, these predictions of the linear sigma model de- 
pend on a single parameter, namely r. The parameter- 
free results in the limit r —> oo (that is, chiral perturba- 
tion theory) are shown in Fig. H]by the thin black lines. 
In fact, adjusting r « 9, one would find that the NJL 
and linear sigma model curves lie exactly on top of each 
other up to \xb about 400 MeV, and remain very close 
in the whole range of chemical potentials displayed in 
Fig. |6l A cautious reader may wonder why this is so, re- 
membering our claim in Sec. lIVI that at /is ~ 229 MeV a 
crossover from the BEC behavior to quark matter occurs. 
The reason for this unexpected agreement even after the 
crossover to quark matter, defined by the change of sign 
of the expression M — is that the quark excita- 

tions will keep a large energy gap. At low temperatures, 
they will not be excited and the bulk thermodynamic 
quantities will be dominated by the condensate (and its 
fluctuations, as we will see later). 

Let us now compare our calculation based on chiral 
modes with the lattice data. A glance at the results pre- 
sented in Ref. [lj| reveals that all three thermodynamic 
observables considered here have substantially different 
values. This is, of course, not so surprising on account 
of the fact that we made no attempt to fit our model 
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parameters to the lattice data. However, this issue de- 
serves more care, and we would now like to explain to 
what extent the lattice data actually can be explained 
with chiral models. 

First of all, one should note that the lattice computa- 
tions are always performed with a nonzero source for the 
diquark operator. While to obtain physical predictions 
one should eventually take the limit of zero source, for the 
sake of comparison it is technically much easier to include 
the nonzero source in the model calculations. Indeed, in 
the NJL model, adding a diquark source amounts simply 
to shifting the diquark field A in the Yukawa interaction 
term in ([5]), and consequently also in the quasifermion 
contribution to the thermodynamic potential (j6|). In the 
linear sigma model Lagrangian (|9|), we would replace Ha 
with Ha+JRe A, J being the diquark source. The result- 
ing gap equations can still be solved analytically, albeit 
implicitly, in the chiral perturbation theory limit, that is, 
t — > oo 2]. Writing the sources H and J in terms of a 
new angle ip as H = J^m^r cos V an d J — f^m^ sin ip, the 
gap equation for 9, which determines the orientation of 
the condensate, reads 



x smv cost* = sm 



•-ip). 



Once this equation is solved for 9, the thermodynamic ob- 
servables are obtained from the generalization of Eq. (fTO)) 
to nonzero source in the limit r — ^ oo, 
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The results of a sample calculation using chiral pertur- 
bation theory with a nonzero diquark source are shown 
in Fig. [7J Since we were not able to fix the input parame- 
ters in the vacuum, we simply fit the prediction of chiral 
perturbation theory to lattice results at nonzero chemi- 
cal potential. We thus obtain ip « 0.10, f v /m v ~ 0.21, 
and run w 0.59 in lattice units. (We used only the data 
for baryon density in the fit, since the pressure may get 
a nontrivial contribution from the gluon sector of QCD. 
Also, only the points with ^is/fn-n J; 2 were taken into 
account.) Doing this, we merely wish to demonstrate 
that a part of the lattice data can be explained using a 
model based just on the chiral symmetry. Obviously, the 
same numerical values can also be obtained with the NJL 
model; we did not do that since it would require tuning 
an even larger number of parameters at will. 

It is clear that the high numerical values of the NJL 
prediction plotted in Fig. [5] do not present a problem. In 
fact, Eq. (fTO)) illustrates convincingly that they scale sen- 
sitively with the ratio j^jm^. What is a problem, how- 
ever, is the large ratio of the heights of the peaks in the 
energy density, and in the pressure and baryon density, 
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FIG. 7. Comparison of lattice results [161 ] (we are indebted to 
S. Hands for providing us with the data) and leading order 
chiral perturbation theory with nonzero diquark source. The 
line notation is the same as in Fig. [S] red solid for pressure 
and green dash-dotted for baryon density. 



as observed on the lattice jTy] . This is the reason why 
we did not plot the energy density in Fig. [7] at all. Both 
the mean-field NJL model and the linear sigma model at 
tree level (and in turn the chiral perturbation theory at 
the leading order) predict all three peaks to have heights 
of roughly the same magnitude. Even though we can- 
not offer and ultimate resolution of this discrepancy, we 
propose a possible explanation. 

Our reasoning goes back to the thermodynamic iden- 
tity e = — p + /iBi^B + Ts. Note that the lattice simu- 
lations were performed at low temperature, roughly an 
order of magnitude smaller than the typical values of 
the chemical potential involved. The SB expressions are 
then well approximated by their zero-temperature limits, 
which satisfy the relation £g B = 3pg B = j^b u sb- The 
SB-normalized thermodynamic quantities (denoted by a 
bar) thus satisfy the identity 



1_ 4_ 
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where the approximation sign refers to the zero- 
temperature limit of the SB values. Using this identity, 
one may from the values of p, e given in [16j determine 
the purely thermal component of the energy density, rep- 
resented by the entropy. One thus finds that Ts/ esb is 
essentially zero for fi B /m v ~ 3 [431 ] . but the sharp peak 
already present in e, of course, remains. We are thus led 
to the conclusion that the high peak in the energy density 
is due to thermal excitations. 

In order to see how this can come about, let us observe 
that the diquark disper sion relat ion below the onset of 
BEC is given by E k = y/k 2 + to 2 -fj, B . Right at the BEC 
transition point, /j-b — fn n , the dispersion is therefore 
quadratic at low momentum, and it still remains rather 
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flat even when the chemical potential increases above m,. 
Moreover, in the whole diquark condensation phase there 
is a massless NG boson stemming from the spontaneous 
breaking of the exact global symmetry associated with 
the conservation of the baryon number. As soon as the 
temperature is nonzero, the NG states will be abundantly 
populated, giving a large entropy, contributing in turn to 
the energy density. Whether this mechanism is physically 
correct can easily be tested by changing the temperature. 
If the peak in the energy density is really mostly due to 
thermal excitations, its height should be very sensitive 
to the changes in temperature. On the other hand, the 
peaks in the pressure and baryon density, not getting a 
contribution from the entropy, should be rather robust, 
being dominated by the condensate contribution. 

With this physical picture in mind, it is now also clear 
how our model calculations should be modified in order 
to take the NG excitations in the account. In the linear 
sigma model or chiral perturbation theory, one must in- 
clude loop effects which introduce the finite temperature 
into the system via the Matsubara summation. Within 
the NJL model, the mean-field approximation only in- 
cludes quasifermionic modes. In order that the NG exci- 
tations be properly taken into account, one must, again, 
go beyond the mean-field approximation. Although both 
these extensions are in principle straightforward (the 
next-to-leading order results within chiral perturbation 
theory being even already available Q), albeit maybe 
somewhat tedious, they may still not yield the desired 
quantitative agreement with the lattice simulations. Pos- 
sible reasons are discussed in the concluding section. 



VII. CONCLUSIONS 

In this paper we investigated the thermodynamics of 
two-color QCD with two light quark flavors. We devised 
a general method to construct interaction Lagrangians of 
the NJL type with the extended global SU(4) symmetry. 
Using the mean-field approximation, we then determined 
the phase diagram in the (fiB,fJ-i) plane, including the 
FF pairing and the isospin splitting of the flavor chiral 
condensates. 

One should emphasize that, since we only considered 
the chiral NJL model, confinement effects are not taken 
into account in our results. This can be done phenomeno- 
logically by adding the Polyakov loop as an additional 
degree of freedom [2(|. However, this will only affect the 
physics at nonzero temperature, since the Polyakov loop 
NJL model reduces to pure NJL at zero temperature. 
Hence, the results of Sec. [TV] will be unaffected. Like- 
wise, Figures [3J HI and [5] will not change qualitatively; 
basically the only difference will be a rescaling of the 
critical temperatures to higher values. 

Finally, in Sec. IVII we discussed the thermodynamics 
of diquark BEC with a particular attention to available 
lattice data. We suggested that a proper inclusion of 
thermal NG excitations is crucial in order to achieve a 



quantitative agreement with them. While it would be 
very nice to pursue our model calculations as far as pos- 
sible, we have some reservations as to whether this would 
lead to the desired goal. The reason is that the lattice 
computations were performed with rather heavy quarks. 
Consequently, the lightest non-NG bosonic state is the 
rho meson with m w /m p ~ 0.8. We see that there is no 
clear separation of scales, and therefore the success of 
any approach based on chiral symmetry and its spon- 
taneous breaking is a priori questionable. One obvious 
remedy to this situation would be to perform lattice sim- 
ulations with lighter quarks, but this may be technically 
rather demanding. On the other hand, one might try to 
avoid using chiral models by resorting to first-principle 
QCD calculations, based, e.g., on the Dyson-Schwinger 
equations 44] or the functional renormalization group 
approach [451 ]. 

Summarizing, we would like to stress once more that 
two-color QCD is a system whose study may teach about 
deconfinement in cold dense matter as well as other as- 
pects of the real world. At present, the techniques that 
would allow us to investigate the phase diagram of two- 
color QCD in detail are on the market. We believe that 
our future efforts will contribute to its general under- 
standing. 
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APPENDIX: THERMODYNAMIC POTENTIAL 
FOR FULDE-FERRELL PAIRING 

Consider for simplicity a nonrelativistic theory of two 
fermion species, ipf,ip^, defined by the (Euclidean) La- 
grangian 

where g is the bare (attractive) coupling and jl the di- 
agonal chemical potential matrix in the flavor space. In- 
troducing a collective field <fr ~ ipiip^, performing the 
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Hubbard-Stratonovich transformation, integrating out 
the fermions, and finally summing over the Matsubara 
frequencies, one arrives at an expression for the mean- 
field thermodynamic potential analogous to Eq. ([5]), 



Q[(t>]= J d 3 x 



Tr{-Jt>(x) + T]og 



(11) 



where 



-4>*(x) 5£ 



is the matrix Hamiltonian in the Nambu space, (ipf, w), 
and the trace operation is to be performed in the operator 
sense. 

So far, no assumption on the particular form of the 
mean field configuration <p(x) was made; in principle this 
is to be determined by minimization of the functional 
51 [4>]. Once again, however, Eq. (fTTj) is only formal be- 
cause TiJ^(x) is badly divergent. Subtracting the same 
expression for the normal phase (that is, <p = 0) at zero 
temperature with the corresponding Hamiltonian J?q, 
one obtains the regularized thermodynamic potential as 



i 



Trj-fjr^-l^ol 



Tlog 1 + e 



-I3.?f?(x) 



(12) 



Within the nonrelativistic theory this prescription can be 
shown to remove all divergences once the coupling g is 
renormalized. (In practice this is usually done by trading 
g for the gap at zero temperature which is physically 
observable and thus free of divergences.) 

Unfortunately, it is impossible to evaluate the trace of 
the log term for a general field configuration <f>{x). The 
FF pairing corresponds to the Ansatz <fi(x) = Ae 2lqx . 
Now the Hamiltonian may already be diagonalized. To 
this end, we perform the unitary transformation, — > 
Jf' = U^JfUq, where U q = diag (e^, e^") and 



Jff" 



f-jV+g) 



-A 



-A 
(-»v-<?) 2 

2m 



The eigenvalues of this Hamiltonian are now easily found 
by Fourier transforming to momentum space. They are 
E k and —E%, where — E k ± 6p, q , 5[i q = Sp — (k ■ 
q)/m, E k = y/t* + A 2 , £ fe = (k 2 + q 2 )/2m - ft, and 
we have defined p = (pf + p\)/2 and Sp = (pf — £t|.)/2. 
Our crucial observation is that since the trace in f| 12|) 
has a well defined mathematical meaning in combination 



with Jtfo, we have to do the same unitary transformation 
U q on it as well. The thermodynamic potential density 
therefore becomes 



n(A,«f) 



V 



A 2 
9 



d 3 k 

(2^)3 

ICfc - SfJ, q \] - 



Ei 



-T^log 



■5tic 



Note that subtracting naively the contribution of the nor- 
mal phase from Eq. (jlip would amount to replacing the 
second term in the curly brackets with | [|^ + Sfx] + \£ k — 
<5/i|] . The difference is finite, but the naive subtraction 
would lead to an extra term in the thermodynamic po- 
tential, independent of A and going to — oo at large q. 
Such a thermodynamic potential would not even have 
a minimum. Our derivation makes sure that the ther- 
modynamic potential is both well defined and physically 
consistent. 

In relativistic systems the same argument can be used, 
but the situation is less simple. The reason is that the 
divergences in the relativistic theory are more severe due 
to different quasiparticle dispersion relation, which is al- 
ready signalized by the presence of a quadratic divergence 
~ A 2 q 2 . The divergences are not completely removed by 
the subtraction of the normal phase and renormalization 
of the coupling constant. This means that the subtrac- 
tion procedure remains somewhat ambiguous and is to 
be understood as a part of the definition of the model. 
This is, however, nothing new in the NJL model where 
one often handles divergent quantities in a rather for- 
mal manner [li|. Our strategy is as follows. We define 
the naive thermodynamic potential, Q(a, p, A(cc)), anal- 
ogously to Eq. (fTTj) . and rewrite it as 

Q(a, p, A(x)) = fifo p, A(x)) - Q(a, p, 0) + Q(a, p, 0). 

The last term is the thermodynamic potential of the 
"normal phase" while in the difference of the first two 
terms the most severe divergences cancel. It is this dif- 
ference that the unitary transformation U q is applied to. 
For the FF pairing this procedure leads to Eq. ([7]). One 
should, however, observe that since the final expression 
is still mildly (logarithmically) divergent, the result de- 
pends on the choice of the term to subtract, which was 
already noted in [l3[ . 

Here we have subtracted the "normal phase" with re- 
spect to the diquark condensate, but with unchanged val- 
ues of the chiral condensates. By dimensional analysis, 
the thermodynamic potential suffers from the following 
g-dependent divergences: a quadratic divergence ~ q 2 
and logarithmic divergences ~ q 2 a 2 ,q 2 p 2 ,q 2 /j, 2 3 j,q 2 A 2 . 
All divergences but the very last one are canceled by our 
subtraction procedure. Fortunately, the remaining di- 
vergence ~ q 2 A 2 is rather mild because in the FF phase 
the gap A is numerically much smaller than in the BCS 
phase. 
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